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INK4 locus at chromosome 9p21 has been reported to be associated with primary open angle glaucoma 
(POAG) and its subtypes along with the associated optic disc parameters across the populations of 
European, Japanese and African ancestries. The locus encodes three tumor suppressor genes namely 
CDKN2A, ARF, CDKN2B and a long non-coding RNA CDKN2B-AS1 (also known as ANRIL). Here, we 
report association study of 34 SNPs from INK4 locus with POAG in a population of Indo-European ancestry 
from the eastern part of India (350 patients and 354 controls). With 81% power to detect genetic association 
we observed only nominal association of rsl011970 (uncorrected p=0.048) with POAG and rsl0120688 
(uncorrected p=0.048) in patients without a high intra-ocular pressure (IOP<21 mm of Hg) compared to 
controls. This study, in contrast to the previous reports, suggests lack of significant genetic association of 
INK4 locus with POAG in East Indian population which needs to be replicated in larger studies in diverse 
world populations. 

Glaucoma is a group of optic neuropathies which is the foremost cause of irreversible blindness. The recent 
estimates from population-based surveys indicate that 60 million people are affected by glaucoma 
worldwide which is estimated to reach 79.6 million by 2020 and India will surpass Europe to have the 
second largest population of glaucoma patients'. In India, the prevalence of blindness is 1% and glaucoma causes 
8.3% of the total blindness^. Primary open angle glaucoma (POAG) is the most common subtype, which is a 
multifactorial complex disease. According to the latest review of the field, family based linkage analyses have 
established 17 loci for POAG out of which only five underlying genes have been identified"". In parallel, case- 
control associations in candidate gene regions as well as genome wide association studies have suggested addi- 
tional loci to be associated with POAG^ 

The INK4 locus at chromosome 9p21 has been found to be associated with POAG and its subtypes along with 
the associated traits like vertical cup-disc ratio (summarised in supplementary Table SI). The association was 
observed in populations of Caucasian, Japanese and some populations of African ancestry. 

The locus encodes three tumor suppressor genes namely CDKN2A, ARF, CDKN2B and a long non-coding 
RNA CDKN2B-AS1 (also known as ANRIL). Recent reviews suggest that CDKN2A and CDKN2B function as cell 
growth regulator controlling cell cycle Gl progression by preventing the activation of cyclin dependent kinases'*. 
It has been reported that ANRIL is involved in the epigenetic regulation of genes at the 9p21 locus'* as well as the 
other genes in trans^. 

Given the apparent importance of this locus in POAG, we tested genetic association of SNPs at this locus in a 
cohort of POAG cases and controls from the state of West Bengal in the eastern part of India. 

Results 

A total of 1392 samples were enrolled for the study. Of these 704 samples (350 cases and 354 controls) were used to 
analyse genotype data for the entire region (Discovery phase) and 688 samples (368 cases and 320 controls) were 
further used to validate candidate tag SNPs from the region. The mean age at onset of the patients recruited in the 
discovery phase was observed to be 48.99 ± 16.6 years and the mean age ofcontrols was 49.75 ± 11.36 years. The 
average lOP of patients was 22.94 ± 7.76 mm Hg. 
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9p21 locus SNPs do not show allelic association in Indian POAG 
patients. Chi-Square analysis did not reveal allelic association of any 
of the 34 SNPs with POAG as represented in regional association plot 
(Fig. 1) except the marginal association of rsl011970. The minor 
allele frequencies and respective p-values are given in Table 1. In 
addition, a marginal association of GG genotype [p = 0.029, 
OR=2.1 (CI= 1.01-4.4)] of rsl412832 was observed as a risk 
towards POAG albeit with a low frequency (0.074 cases/0.036 
controls). Our study had >80% power to detect an allelic 
association for this region. The details of the power analysis along 
with comparison with other studies reporting positive association is 
summarized in Table 2. To further confirm the negative association 
of 9p21 locus in Indian POAG patients we selected three tag SNPs 
from this region and genotyped them in an independent cohort of 
368 cases and 320 controls from the same population background. 
We did not observe any allelic association (Table 3). To rule out any 
possible error in the analysis we also analysed publically available 
GLAUGEN data for this locus which reported positive association 
and could reproduce association in 12 SNPs from this locus (data not 
shown). 

Genotype association analysis with cup-disc ratio. Among the 
markers examined at this locus, rs 1063 192 has been shown to be 
associated with reduced cup-disc ratio in Caucasian and Japanese 
population backgrounds (supplementary Table SI). We tested 
association of cup-disc ratio with the alleles and genotypes of this 
SNP in 215 POAG cases by linear regression analysis but did not find 
any association (data not shown). 

Selection of tag SNPs and haplotype association. Further, for 
haplotype analysis of this locus, 20 SNPs were chosen as tag SNPs 
on the basis of their pairwise r^ values in our dataset (r^>0.8) 
(Supplementary Fig. 1). Thirteen haplotypes were generated 



among which only one showed significant frequency difference 
between cases and controls although the overall frequency of 
occurrence was low in both groups (0.005 cases/0.017 controls) 
(Supplementary Table S2). 

Association analysis in HTG and Non-HTG subjects. Many studies 
have reported association of this locus with POAG patients who do 
not suffer from high lOP (Supplementary Table S 1 ) . In accordance to 
these studies, we divided our patients into two subgroups: (a) high 
tension glaucoma (HTG) cases having IOP>21 mm of Hg in at least 
one eye (n = 232), and (b) non-HTG cases with IOP<21 mm of Hg 
in both eyes (n=104) and compared the allele frequencies of SNPs 
separately against 354 controls. The average lOP estimate of HTG 
subgroup was 25.63 ± 7.79 mm Hg and of non-HTG subgroup was 
16.98 ± 2.71 mm Hg. In non-HTG subgroup analysis, rsl0120688 
was observed to have a marginal association (p<0.048); but none 
were found to be associated with the HTG subgroup (Table 1). 

Discussion 

Our study suggests a negative association of the INK4 locus with 
POAG in Indian sub-continent. There are very marginal/weak 
associations of (a) rsl011970 with POAG, and (b) rsl0120688 with 
non-HTG subgroup of the disease which will not sustain statistical 
correction for multiple testing. The association of rslO 120688 with 
normal-tension glaucoma was found in the meta-analysis of 
GLAUGEN and NEIGHBOR study populations by Wiggs et al 
(Supplementary table SI) and also in African-Americans, however 
it did not pass multiple-testing correction*". Burdon et al also found 
association of this SNP with lOP (Supplementary table SI). The 
genetic association might be found positive if a larger study is under- 
taken with exclusively normal tension glaucoma patients. 

It is worthwhile to mention here that the 9p21 locus has been 
demonstrated to have remarkable variability in its genetic association 
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Figure 1 | Regional association plot of INK4 locus in East Indian population. This plot shows association results of studied SNPs from 9p21 region with 
POAG in East Indian cohort. The circles represent the genotyped SNPs. The horizontal axis is the genomic context of the region studied (along with the 
genes). The left vertical axis represents negative logarithm of p-value and the right vertical axis is the recombination frequency of the region. Purple circle 
(marked by arrow) represents rslOl 1970 which has minimum p-value. Overall the data suggests minimal association (if any) of this locus 
with POAG. 
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Table 1 | Details of INK4 locus SNPs at 9p21 in East Indian cohort 

Controls 



SNP 



r$2811708 

rs3731239 

rs2811709 

rs4074785 

rs3731217 

rs32 18020* 

rs2811712 

rs3218018 

rs32 17992 

rsW63192* 

rs2069426 

rs32 17980 

rs2069422 

rs573687* 

rsl 3298881 

rs545226 

rs7049105* 

rs2157719* 

rsl 75941 7 

rs2151280* 

rsl412829* 

rsl 333034 

rsl 01 20688* 

rsl011970 

rs4977756* 

rs7855162 

rsl 41 2832* 

rsl Oil 6277* 

rsl 547705 

rsl 333040 

rs4977574 

rsl 0965235 

rs944797 

rs2383207* 



Minor 
allele _ 

A 
G 
A 
A 
C 
A 
G 
C 
A 
G 
A 
A 
C 
A 
G 
G 
A 
G 
A 
G 
G 
G 
G 
A 
G 
G 
G 
C 
C 
G 
A 
A 
A 
A 



(n=708j POAGln=700) 
^AAF ^AAF p-value OR(95%CI) ^AAF 



HTG** (n=464l 
p-value OR(95%CI) 



Non-HTG** ln=208) 
MAF p-value OR(95%CI) 



0.17 

0.23 

0.05 

0.17 

0.1 1 

0.37 

0.15 

0.06 

0.39 

0.25 

0.06 

0.01 

0.06 

0.23 

0.02 

0.39 

0.4 

0.25 

0.03 

0.4 

0.24 

0.06 

0.39 

0.27 

0.24 

0.09 

0.22 

0.38 

0.03 

0.36 

0.48 

0.08 

0.46 

0.38 



0.17 
0.22 
0.04 
0.19 
0.13 
0.34 
0.13 
0.04 
0.42 
0.24 
0.04 
0 

0.04 
0.23 
0.02 
0.42 
0.37 
0.24 
0.02 
0.38 
0.23 
0.04 
0.35 
0.23 
0.24 
0.09 
0.24 
0.36 
0.03 
0.35 
0.45 
0.07 
0.43 
0.36 



0.96 

0.79 

0.38 

0.32 

0.46 

0.07 

0.23 

0.16 

0.21 

0.9 

0.09 

0.13 

0.12 

0.86 

0.83 

0.31 

0.35 

0.76 

0.28 

0.39 

0.85 

0.15 

0.07 

0.048 

0.95 

0.76 

0.38 

0.51 

0.56 

0.77 

0.16 

0.41 

0.23 

0.44 



0.99(0. 
0.97(0. 
0.79(0. 
1.15(0. 
1.13(0. 
1.22(0. 
0.83(0. 

0.7(0. 
1.15(0. 
0.98(0. 
0.66(0. 
0.38(0. 
0.69(0. 
1.02(0. 
1.08(0. 
1.12(0. 

0.9(0. 
0.96(0. 
0.67(0. 
0.91(0. 
0.98(0. 

0.7(0. 
0.82(0. 
0.78(0. 
0.99(0. 
0.94(0. 
1.12(0. 
0.93(0. 
0.83(0. 
0.97(0. 
0.86(0. 
0.85(0. 
0.88(0. 
0.92(0. 



75-1.31) 
75-1 .24) 
47-1.33) 
88-1.51) 
82-1.55) 
98-1.51) 
62-1.12) 
43-1.15) 
93-1 .42) 

77- 1 .26) 
4-1 .08) 
1-1.43) 

42- 1.1 1) 

8- 1.31) 
53-2.2) 

9- 1.38) 
73-1.12) 

75- 1.23) 
32-1.4) 

73- 1.13) 

76- 1 .25) 

43- 1.15) 
66-1 .02) 
62-1) 

78- 1.27) 
65-1.36) 
87-1 .43) 
75-1.15) 
45-1 .54) 
78-1.21) 
7-1 .06) 
58-1.25) 
71-1.09) 

74- 1.14) 



0.16 

0.22 

0.034 

0.18 

0.12 

0.41 

0.13 

0.038 

0.42 

0.24 

0.036 

0.004 

0.038 

0.23 

0.02 

0.42 

0.37 

0.24 

0.017 

0.37 

0.23 

0.038 

0.35 

0.23 

0.24 

0.086 

0.24 

0.37 

0.028 

0.35 

0.45 

0.073 

0.44 

0.36 



0.90 
0.99 
0.31 
0.43 
0.66 
0.1 1 
0.24 
0.14 
0.28 
0.94 
0.08 
0.20 
0.12 
0.79 
0.97 
0.36 
0.41 
0.84 
0.35 
0.48 
0.93 
0.14 
0.20 
0.1 1 
0.72 
0.81 
0.50 
0.81 
0.66 
0.95 
0.34 
0.48 
0.47 
0.73 



0.98(0. 

1(0. 
0.73(0. 
1.13(0. 
1 .08(0. 
1.21(0. 
0.82(0. 
0.66(0. 
1.14(0. 
1.01(0. 

0.6(0. 
0.38(0. 
0.64(0. 
1 .04(0. 
1 .02(0. 
1.12(0. 

0.9(0. 
0.97(0. 
0.67(0. 
0.92(0. 
0.99(0. 
0.66(0. 
0.85(0. 

0.8(0. 
1 .05(0. 
0.95(0. 

1 . 1 (0. 
0.97(0. 
0.86(0. 
0.99(0. 
0.89(0. 
0.85(0. 
0.92(0. 
0.96(0. 



72-1 .34) 
76-1 .33) 
4-1 .34) 
83-1.53) 

76- 1 .56) 
95-1 .54) 
58-1.15) 
37-1.16) 
9-1 .45) 

77- 1.33) 
34-1 .07) 
08-1 .79) 
36-1.13) 
79-1 .37) 
45-2.29) 
88-1.42) 

71- 1.15) 

74- 1 .28) 
29-1 .56) 

72- 1.17) 

75- 1 .3) 



37-1 
67-1 
61-1 



16) 
09) 
05) 



8-1.38) 
63-1 .44) 
83-1.45) 
76-1 .24) 
43-1.71) 
78-1.27) 

71- 1.13) 
55-1.32) 

72- 1.16) 
75-1.22) 



0.16 

0.19 

0.028 

0.19 

0.13 

0.43 

0.12 

0.033 

0.44 

0.22 

0.033 

0.004 

0.033 

0.22 

0.028 

0.43 

0.35 

0.23 

0.019 

0.35 

0.22 

0.033 

0.31 

0.21 

0.21 

0.086 

0.24 

0.33 

0.028 

0.33 

0.42 

0.067 

0.41 

0.34 



0.93 

0.38 

0.26 

0.45 

0.33 

0.08 

0.3 

0.17 

0.16 

0.56 

0.15 

0.4 

0.15 

0.78 

0.52 

0.25 

0.28 

0.6 

0.6 

0.25 

0.63 

0.17 

0.048 

0.08 

0.41 

0.86 

0.63 

0.25 

0.79 

0.5 

0.16 

0.42 

0.2 

0.4 



0.98(0, 
0.84(0, 
0.61(0, 
1.17(0, 
1.25(0, 
1.32(0, 
0.79(0, 
0.57(0, 
1.25(0, 

0.9(0, 
0.55(0, 
0.42(0, 
0.55(0, 
0.95(0, 
1.37(0, 

1.2(0, 
0.84(0, 
0.91(0, 
0.75(0, 
0.83(0, 
0.91(0, 
0.57(0, 
0.72(0, 
0.72(0, 
0.86(0, 
0.95(0, 
1.09(0, 
0.83(0, 
0.88(0, 
0.89(0, 

0.8(0, 
0.78(0, 
0.82(0, 
0.87(0, 



65- 1.48) 
57-1.24) 
25-1.47) 
78-1.74) 
8-1.98) 
97-1.81) 
5-1 .24) 
25-1.28) 
91-1.71) 

62- 1 .29) 
24-1.25) 

05- 3.4) 

24- 1.25) 

66- 1.38) 
53-3.58) 
88-1 .64) 
61-1.16) 

63- 1.31) 

25- 2.23) 

6- 1.14) 

63- 1.32) 
25-1.28) 
52-1) 
49-1.04) 
59-1.25) 
55-1.65) 
76-1.57) 
6-1.15) 
36-2.2) 

64- 1.24) 
59-1 .09) 
43-1.43) 
6-1.12) 
63-1.2) 



*SNPs which have been found to be asBocialed with POAG in the previous studies, 
**HTG and Non-HTG cohorts are the subtypes of POAG cohort. 

Note: p-values in bold represents marginal association of SNPs. Italicized rs-IDs represent 20 tag SNPs chosen for haplotype association test shown in supplementary Table 52; Control data is same for POAG, 
htTG and Non-hlTG groups. 

Abbreviations- n: Number of chromosomes, htTG: High Tension Glaucoma, Non-HTG: Non-High Tension Glaucoma, MAF: Minor allele frequency, OR: Odds ratio, O: Confidence interval. 



and potential role in POAG pathogenesis. For example, some GWA 
studies identified this locus to be significantly associated 
(Supplementary Table SI), while other studies carried on similar 
population did not support the association'. Further, a study on 
Afro-Caribbean ancestry claimed association of this locus" which is 
supported by similar observation made in African- Americans 
though with marginal association*". However, in the same study the 
association did not sustain when investigated in populations from 
Ghana, West Africa*". Similarly, a study in Brazilian population on 
one SNP (rs4977756) from this locus also reported lack of any asso- 
ciation with POAG". 

It is possible that the effect of this locus on the disease might be 
exerted by, yet unknown, genetic or environmental factor(s) that are 



variable in different population groups and not being considered 
in any of the current studies undertaken to address POAG 
pathogenesis. 

Recently, Janssen et al found that among the reported 65 POAG 
candidate genes, which includes 5 familial genes, 33 highly likely 
candidate genes and 27 less-likely candidate genes for POAG belong 
to only four functional molecular networks'. Genes at 9p21 locus are 
part of their molecular network '4' constructed by Ingenuity know- 
ledge database, which consisted of 12 candidate POAG genes where 
main functionalities were ophthalmic disease, cell cycle, connective 
tissue development and function. It might be possible that in the 
populations where this locus has not been found to be associated, 
other molecular networks play major roles but it will require more 



Table 2 | Power calculations 


to detect Genetic association 






SNP Allele 


Wiggs etai. 2012 
(31 46 cases, 3487 controls) 


Fan B.J. etal2011 
(536 cases, 336 controls) 


East Indian population 
(350 cases, 354 controls) 




POWER (%} Allele frequency 


POWER {%} Allele frequency 


POWER (%l Allele frequency 


rsl063192 A 


100 0.57 


90.32 0.58 


81.21 0.75 


Note: Power calculations were done for 'A' 
frequency were taken from respective slue 


allele by keeping relative risk as 1.8 ( Homozygous genotype) and 1 . 2 (Heterozygous genotype) calculated from genotype data given hy Fan BJ.et al, 20 1 1 and allele 
es (Supplementary Table S 1 ). 
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Table 3 | Details of association of three replicated tag SNPs representing three different LD blocks 

Discovery phase Replication phase Combined 

, 350 cases, 354 controls 368 cases, 320 controls 71 8 cases, 674 controls 

SNP Minor allele 

MAF(Cases/Controlsj p-value MAF(Cases/Confrols) p-vaiue MAF(Cases/Contro!sj p-va!ue 

rs32 1 7992 A 0.42/0.39 0.20 0.39/0.4 0.74 0.40/0.39 0.5 1 

rs4977756 G 0.24/0.24 0.94 0.27/0.25 0.38 0.25/0.24 0.53 

rslOl 16277 C 0.36/0.38 0.50 0.41/0.39 0.39 0.39/0.38 0.84 



extensive genetic analysis supported by gene expression assays and 
functional characterization. 

Our finding of lack of genetic association of INK4 locus in a 
modest cohort of East Indian POAG patients indicates biological 
complexity of this region. This requires further validation by under- 
taking a much larger study population; but the underlying variable 
observations made in different populations across the globe remains 
a greater challenge to resolve. 

Methods 

Selection criteria of study individuals. Diagnosis of patients has been done in Dristi 
Pradip eye clinic, Jodhpur Park, Kolkata. The patients have been diagnosed through 
clinical ocular and systemic examinations. The inclusion and exclusion criteria for 
samples were the same as reported earlier*". Briefly, the patients were recruited if they 
were positive for 2 out of the 3 criteria, namely, Intra-ocular pressure (IOP)>21 mm 
of Hg, glaucomatous field damage and significant cupping of the optic disc. 
Individuals with ocular hypertension and with any history of inflammation or ocular 
trauma (past & present) were excluded from this study. 

Controls were selected from the general population without any history of ocular 
disease and who tested negative for POAG by means of routine eye examination for 
glaucoma as described above. The study protocol adhered to the tenets of the 
Declaration of Helsinki and was approved by the Institutional Review Board. 

Collection of blood samples and genomic DNA preparation. Peripheral blood was 
collected with EDTA from the POAG patients and the individuals without POAG 
with their written consent. Genomic DNA was prepared from fresh whole blood 
using the PAX gene blood DNA isolation kit (Qiagen, HUden, Germany) according to 
the manufacturer'sprotocoI.TheDNA was dissolved in TE{10 mMTris-HCl, 1 mM 
EDTA, pH 8.0). 

Genotyping. The genotype data of 38 SNPs were extracted from a genome-wide 
genotyping dataset for 350 POAG cases and 354 controls which was generated using 
lUumina Human660W-Quad chip {Illumina, Inc., San Diego, CA, USA). Four SNPs 
were removed in QC for minor allele frequency (MAF<0.01) and Hardy- Weinberg 
equilibrium {HWE p<0.01). Chi-square association was performed on 34 SNPs 
spanning 142.53 kb in reference human genome (hgl8). The SNP and the sample call 
rates were >95% and >98%, respectively. The selected tag SNPs for validation were 
genotyped using the GoldenGate genotyping assay (lUumina GGGT assay). 

Statistical analysis. Power calculations were done for 'A' allele of rsl063192 which is 
the most replicated SNP in different reported studies by Genetic power calculator 
software (http://pngu.mgh.harvard.edu/~purcell/gpc/)'\The relative risk was kept 
constant as 1.8 (for homozygous genotype) and 1.2 (for heterozygous genotype) as 
reported hy Fan B.J. etal (Supplementary Table SI) and allele frequency of respective 
populations were taken as parameters. Chi-square test of association was performed 
using PLINK software*^. The regional plots were generated by an online available 
Locus Zoom tool (http://csg.sph.umich.edu/locuszoom/) LD calculations and 
haplotype association were performed using haploview software version 4^*. 
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